CLNS 10/2068 
CLEO 10-05 

Search for rare and forbidden decays of charm and 
charmed- strange mesons to final states h ± e^e + 

P. Rubin, 1 N. Lowrey, 2 S. Mehrabyan, 2 M. Selen, 2 J. Wiss, 2 J. Libby, 3 M. Kornicer, 4 

R. E. Mitchell, 4 M. R. Shepherd, 4 C. M. Tarbert, 4 D. Besson, 5 T. K. Pedlar, 6 
J. Xavier, 6 D. Cronin-Hennessy, 7 J. Hietala, 7 P. Zweber, 7 S. Dobbs, 8 Z. Metreveli, 8 
K. K. Seth, 8 A. Tomaradze, 8 T. Xiao, 8 S. Brisbane, 9 L. Martin, 9 A. Powell, 9 
P. Spradlin, 9 G. Wilkinson, 9 H. Mendez, 10 J. Y. Ge, 11 D. H. Miller, 11 I. P. J. Shipsey, 11 
B. Xin, 11 G. S. Adams, 12 D. Hu, 12 B. Moziak, 12 J. Napolitano, 12 K. M. Ecklund, 13 

J. Insler, 14 H. Muramatsu, 14 C. S. Park, 14 L. J. Pearson, 14 E. H. Thorndike, 14 
F. Yang, 14 S. Ricciardi, 15 C. Thomas, 9 ' 15 M. Artuso, 16 S. Blusk, 16 R. Mountain, 16 
T. Skwarnicki, 16 S. Stone, 16 J. C. Wang, 16 L. M. Zhang, 16 G. Bonvicini, 17 D. Cinabro, 17 
A. Lincoln, 17 M. J. Smith, 17 P. Zhou, 17 J. Zhu, 17 P. Naik, 18 J. Rademacker, 18 
D. M. Asner, 19 '* K. W. Edwards, 19 K. Randrianarivony, 19 G. Tatishvili, 19 - * R. A. Briere, 20 
H. Vogel, 20 P. U. E. Onyisi, 21 J. L. Rosner, 21 J. P. Alexander, 22 D. G. Cassel, 22 
S. Das, 22 R. Ehrlich, 22 L. Fields, 22 L. Gibbons, 22 S. W. Gray, 22 D. L. Hartill, 22 
B. K. Heltsley, 22 D. L. Kreinick, 22 V. E. Kuznetsov, 22 J. R. Patterson, 22 D. Peterson, 22 
V ■ D. Riley, 22 A. Ryd, 22 A. J. Sadoff, 22 X. Shi, 22 W. M. Sun, 22 and J. Yelton 23 

j§ ; (CLEO Collaboration) 

1 George Mason University, Fairfax, Virginia 22030, USA 
2 University of Illinois, Urbana- Champaign, Illinois 61801, USA 
3 Indian Institute of Technology Madras, 



O 

m 

(N 



> 



O ■ Chennai, Tamil Nadu 600036, India 

4 Indiana University, Bloomington, Indiana 47405, USA 
5 University of Kansas, Lawrence, Kansas 66045, USA 
q ■ 6 Luther College, Decorah, Iowa 52101, USA 

7 University of Minnesota, Minneapolis, Minnesota 55455, USA 
8 Northwestern University, Evanston, Illinois 60208, USA 
9 University of Oxford, Oxford 0X1 3RH, United Kingdom 
^ ■ 10 University of Puerto Rico, Mayaguez, Puerto Rico 00681 

11 Purdue University, West Lafayette, Indiana 47907, USA 
12 Rensselaer Polytechnic Institute, Troy, New York 12180, USA 
13 Rice University, Houston, Texas 77005, USA 
14 University of Rochester, Rochester, New York 14627, USA 
15 STFC Rutherford Appleton Laboratory, Chilton, 
Didcot, Oxfordshire, 0X11 OQX, United Kingdom 
16 Syracuse University, Syracuse, New York 13244, USA 

17 Wayne State University, Detroit, Michigan 48202, USA 

18 University of Bristol, Bristol BS8 1TL, United Kingdom 
19 ' Carleton University, Ottawa, Ontario, Canada K1S 5B6 

20 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
21 University of Chicago, Chicago, Illinois 60637, USA 
22 Cornell University, Ithaca, New York 14853, USA 
23 University of Florida, Gainesville, Florida 32611, USA 
(Dated: September 9, 2010) 



1 



Abstract 

We have searched for flavor-changing neutral current decays and lepton-number-violating decays 
of D + and Df mesons to final states of the form /i =t e T e + , where h is either 7r or K. We use 
the complete samples of CLEO-c open-charm data, corresponding to integrated luminosities of 818 
pb _1 at the center-of-mass energy Eqm = 3.774 GeV containing 2.4 x 10 6 D + D~ pairs and 602 pb _1 
at E CM = 4.170 GeV containing 0.6 x 10 6 D^Df pairs. No signal is observed in any channel, and 
we obtain 90% confidence level upper limits on branching fractions B{D + — > ir + e + e~) < 5.9 x 10~ 6 , 
B{D+ -> ir-e+e + ) < l.lxlO" 6 , B{D+ -»• K + e + e~) < 3.0xl0~ 6 , B(D+ -> K"e+e + ) < 3.5xl0- 6 , 
B(D+ -> vr+e+e") < 2.2 x 10" 5 , B{D+ -> i:~e + e + ) < 1.8 x 10" 5 , B(D+ -> K + e + e~) < 5.2 x 10~ 5 , 
and B(D+ -> ^- e +e+) < 1.7 x 10~ 5 . 



* Present address: Pacific Northwest National Laboratory, Richland, WA 99352 
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I. INTRODUCTION 



As an extension of our previously reported [1] search for rare and forbidden decays of the 
D + charm meson, D + — > h ± e T e + , we report an analysis using CLEO-c's full open-charm 
data sample for D + , and also a search for — > h ± e T e + with CLEO-c's full data 
sample. Here, h is either it or K, and charge-conjugate modes are implicit throughout this 
article. These decays probe flavor-changing neutral currents (FCNC), in D + — > 7i + e + e~ and 
D+ — > K + e + e~, and lepton number violations (LNV), in D + — > h~e + e + and D+ — > h~e + e + . 
These decays are either highly suppressed or forbidden in the standard model (SM), but can 
be significantly enhanced by some non-SM physics scenarios [2-7]. Standard model short- 
distance FCNC decays are expected to be of order 10~ 10 to 10~ 9 [3, 5], but long-distance 
vector-pole induced decays of D + or £)+ — > h + V° — > h + e + e~ (where V° is an intermediate 
vector meson p°, u, or 0) are expected to be of order 1CT 6 to 1CT 5 [3, 5]. To observe an 
enhancement in FCNC due to non-SM physics, we need to search for dielectron mass regions 
away from the vector poles. Measuring long-distance induced decay itself might be helpful 
to understand the long-distance dynamics in the b sector, such as inclusive b — > 57 decay or 
exclusive B — > p*y and B — > K*^ decays related to extracting Cabibbo-Kobayashi-Maskawa 
matrix elements |V^( djS )|. On the other hand, observation of LNV (AL = 2) decays could be 
an indication of a Majorana nature of neutrinos [6, 7]. 

We have used two sets of open-charm data samples collected by the CLEO-c detector 
in e + e~ collisions provided by the Cornell Electron Storage Ring (CESR). The integrated 
luminosities are 818 pb" 1 at the center-of-mass energy E C m = 3.774 GeV near the peak of 
the ■0(3770) resonance which decays to DD pairs, and 602 pb -1 at Ecm = 4.170 GeV near 
the peak of D^DJ pair production. The 3.774 GeV data set contains 2.4 x 10 6 D + D~ 
pairs and is used to study D + — >■ h ± e T e + decays. The 4.170 GeV data set contains 0.6 x 10 6 
D^Df pairs, and is used to study D+ — > h ± e T e + decays. 

The remainder of this article is organized as follows. The CLEO-c detector is described in 
Sec. II. Event selection criteria are described in Sec. III. Features of background processes, 
our suppression strategy, and signal sensitivity are discussed in Sec. IV. Results are pre- 
sented as plots and tables in Sec. V. Systematic uncertainties associated with the branching 
fractions and their upper limits are discussed in Sec. VI. Finally, a summary of our results 
with systematic uncertainties is provided in Sec. VII. 

II. THE CLEO-c DETECTOR 

The CLEO-c detector [8-11] is a general-purpose solenoidal detector equipped with four 
concentric components: a six-layer vertex drift chamber, a 47-layer main drift chamber, a 
ring-imaging Cherenkov (RICH) detector, and a cesium iodide electromagnetic calorimeter, 
all operating inside a 1 Tesla magnetic field provided by a superconducting solenoidal mag- 
net. The detector provides acceptance of 93% of the full 4-7T solid angle for both charged 
particles and photons. The main drift chamber provides specific- ionization (dE/dx) mea- 
surements that discriminate between charged pions and kaons. The RICH detector covers 
approximately 80% of 4n and provides additional separation of pions and kaons at mo- 
mentum above 700 MeV. Hadron identification efficiencies are approximately 95% with 
misidentification rates of a few percent [12]. Electron identification is based on a likelihood 
variable that combines the information from the RICH detector, dE/dx, and the ratio of 
electromagnetic shower energy to track momentum (E/p). Typical electron identification 
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efficiency is well over 90% on average with the pion fake rate less than 0.1% and the kaon 
fake rate less than a percent [13, 14]. 

A GEANT-based [15] Monte Carlo (MC) simulation is used to study efficiencies of signal 
and background events. Physics events are generated by evtgen [16], tuned with improved 
knowledge of charm decays, and final-state radiation (FSR) is modeled by photos [17]. 
Nonresonant FCNC and LNV signal events are generated according to phase space. 

III. EVENT SELECTION 

Signal candidates are formed from sets of well-measured drift chamber tracks consistent 
with coming from the nominal interaction point. Charged pions and kaons are identified 
from the tracks with momentum greater than 50 MeV and with | cos 6*| < 0.93, where 9 
is the angle between the track and the beam axis. Electron candidates are required to be 
above 200 MeV with | cos#| < 0.90 to ensure that E/p is well measured. 

At E C m = 3.774 GeV, for each signal candidate of the form D + — > h ± e T e + (where 
h is either n or K), two kinematic variables are computed to define a signal region: the 
energy difference AE = E^+ — -Ebcam and the beam-constrained mass difference AM^ C = 
[-^bcam - P 2 d+] 1 ^ 2 ~ m D+i where (Ed+, Pd+ ) is the four- momentum of the signal D + candidate, 
-Kbeam is the beam energy, and m D + is the nominal [18] mass of the D + meson. To improve 
the resolution of the kinematic variables, we recover bremsstrahlung photon showers within 
100 mrad of the direction of the electron candidates. We define a signal box for further 
analysis as (AE, AM^ C ) = (±20 MeV, ±5 MeV), which corresponds to about 3-standard 
deviations of the kinematic variables. Because the expected contribution from the resonant 
decay B(D + — > (pn + — > n + e + e~) ~ (9(10 -6 ) is within our sensitivity, we further subdivide 
D + — > 7r + e + e~ candidates into two channels: resonant D + — > <p(e + e^)n + and nonresonant 
D + — > 7r + e + e~ for the FCNC search. If the dielectron invariant mass M ee of the signal 
candidate is within ±20 MeV of the nominal [18] mass of the meson, we treat it as a 
resonant D + — > (f)(e + e~)n + candidate and exclude it from the D + — >■ 7r + e + e~ candidates. 

Similarly, at E C m = 4.170 GeV, for each signal candidate of the form Df — > h ± e T e + , the 
following two variables are computed to define a signal region: the mass difference AM = 
M D + — m D + and the recoil mass (against the signal candidate) difference AM rccoi i(D^) = 

[(E — E D +) 2 — (p — Pd+) 2 ] 1 ^ 2 — m D*+' where M D + is the invariant mass of the signal 
candidate, m D + is the nominal [18] mass of the D+, (E , p ) is the total four-momentum of 
the e + e~ beam taking the finite beam crossing angle into account, (E D +,p D +) is the four- 
momentum of the signal candidate with E D + = [m^ + ±p^ + ] 1 / 2 , and m D *+ is the nominal [18] 

s J_J S L) s s 

mass of the D* + . The same bremsstrahlung recovery is performed and the — > 7i + e + e~ 
channel is subdivided into resonant 0(e + e~)7r + and nonresonant channels. The signal box 
is defined as (AM, AM rccoil (£>+)) = (±20 MeV, ±55 MeV) for further analysis. The broad 
recoil mass window ±55 MeV is required to allow both primary and secondary (from D* + — > 
£>+7 or D* s + ->■ -D+7T ) £>+ candidates to be selected. 

IV. ANALYSIS 

Backgrounds are dominantly from events with real electrons, particularly from D semilep- 
tonic decays. The majority of combinatorial background events are from double charm 
semileptonic decays, typically 4 or less charged particles in the event with large missing 
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energy due to the missing neutrinos. Hadronic decays involving 7-conversion and 7r° (77, oj) 
Dalitz decay, or accompanied by another charm semileptonic decay, can mimic the h ± e T e + 
signal, as well. Because of the low probability of hadrons being misidentified as electrons [13], 
background from DD decays to 3-body charged-particle hadronic decays (such as K~it + -n + , 
7r~7r + 7r + , K S K + , K + K~n + ) are negligible after two electrons are identified, and they do not 
peak at the signal region due to the wrong mass assignments for the hadrons misidentified as 
electrons. That is, DD backgrounds are predominantly associated with the semileptonic de- 
cays and non-DD (qq continuum, r-pair, radiative return, or QED events) backgrounds are 
associated with the 7-conversion and Dalitz decays. All of these backgrounds are nonpeaking 
or peak away from the signal regions. 

Our background suppression criteria tuning procedure for D + — > h ± e T e + channels is 
detailed in our previous article [1]. We have used the same background rejection criteria 
with the four kinematic variables to reject the above-mentioned backgrounds in D + channels 
and revised the criteria to accommodate the Df channels. The other side total energy -Esther 
is the sum of energies of all particles other than those making up the signal candidate. We 
use this variable to reject events associated with semileptonic decays, mainly for double 
charm semileptonic decays, in which the visible other side energy would be small due to the 
undetectable missing neutrinos. We reject candidates if -Eothcr < 1-0 GeV for D + — > 7r + e + e~, 
£ thcr < 1-3 GeV for D + K + e + e~, E othcv < 1.4 GeV for D+ -> vr+e+e", and £ other < 
1.7 GeV for Df — > K + e + e~ . For the LNV modes, we reject candidates if the number 
of tracks in the event is 4 or fewer and -E thcr < 0.5 GeV. Semileptonic events involving 
K s — > 7T + 7r - in the final state can mimic the signal in 7r + e + e~ channels. We have used the 
invariant mass M n + n - to veto these events. We veto the candidate when the charged pion 
in the signal candidate combined with any other unused oppositely charged track satisfies 
\M n + n - — m K o \ < 5 MeV, where m K o is the nominal [18] mass of the Kg. Real electrons 
from 7-conversion and Dalitz decays are suppressed by using the dielectron invariant mass 
squared q 2 computed from the signal electron positron pair, or ?o ther computed using one 
signal side electron (positron) combined with any oppositely charged unused track. We 
veto candidates if q 2 < 0.01 GeV 2 or g 2 ther < 0.0025 GeV 2 . For Df, we have required 
the solo photon from D* + decays to Df^f to be explicitly reconstructed to further suppress 
underlying nonstrange-charmed meson backgrounds at E C m = 4.170 GeV, by requiring the 
recoil mass of the signal candidate plus solo photon M rccoi i(.D+ + 7) to be within ±30 MeV of 
the nominal [18] Df mass. Regardless of whether the signal candidate is the primary or 
secondary Df, for the decay e + e~ — > D^DJ — > (Df^y)Df, the mass of the system recoiling 
against the D+ plus 7 should peak at the D+ mass. 

The analysis was done in a blind fashion. Before we opened the signal box, all above- 
mentioned criteria were optimized using MC events with a sensitivity variable which is 
defined as the average upper limit one would get from an ensemble of experiments with the 
expected background and no signal, 

X]jV obs =0 ^(^obs| V cxp )P(iV b s |iVexp) 

S ~ We ' (1) 

where iV ex p is the expected number of background events, iV bs is the observed number 
of events, C is the 90% confidence coefficient upper limit on the signal, V is the Poisson 
probability, iV is the number of D + or Df, and e is the signal efficiency In addition to 
the signal MC samples, four types of background MC samples are utilized to optimize the 
background suppression criteria: 20 times the data sample for open-charm (DD, D*D, 
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FIG. 1. Scatterplots of AM^ C vs AE. The two contours for each mode enclose regions determined 
with signal MC simulation to contain 50% and 85% of signal events, respectively. The signal region, 
defined by (AE, AM bc ) = (±20 MeV, ±5 MeV), is shown as a box. 



D*D*, D*Dir D^D~, and D^ + D~), 5 times the data sample of noncharm uds continuum 
(qq), r-pair, and radiative return to the ip(2S). To normalize background MC events to 
match the expected number of the data events, we have used integrated luminosity and 
cross sections for each process. For D + — > h ± e T e + events at Eqm = 3774 MeV, we have 
used (Jd+d- = 2.91 nb [12], ct d oqo = 3.66 nb [12], a q q = 13.9 nb [19], cr r + T - = 3.0 nb 1 , and 
radiative return to the ip(2S) <trr = 3.4 nb [20]. For Df — > h ± e T e + events at Eqm = 4170 
MeV, we have used a D *± D T = 0.916 nb [21] (and used other open-charm cross sections from 
the same reference), a q q = 11.4 nb [19], <J T + T - = 3.6 nb, and radiative return to the ip(2S) 
crr — 0.50 nb [20]. We have found that the agreements between data and MC simulated 
events are excellent in various kinematic variables used in the background suppression, giving 
us confidence in our optimization procedure using our MC samples. Possible systematic 
uncertainties due to the data and MC differences are assessed in Sec. VI. 



V. RESULTS 



Scatterplots of AE vs AMt, c and AM(D^) vs AM reco ii(D+) for signal candidates with 
all background suppressions applied are shown in Figs. 1 and 2. Except for the 0(e + e~)7r + 
channels, we find no evidence of signals, and we calculate 90% confidence level upper lim- 
its (UL) on the branching fractions based on Poisson processes with background [22] (e.g. 
Section 28.6.4 Poisson processes with background therein) as summarized in Table I: 

UL = »J. (2) 

For D + and — > <j)(e + e~)ir + channels, we find weak evidence of signals with significance 
3.5 for the D + and 1.8 for the D+, so both branching fractions and upper limits are shown 
in Table I. 



1 With the lowest-order QED calculation, cr(e + e~ -> t+t~) = 27ra 2 /3(3 - /3 2 )/(3s), where = (1 - 
4m 2 ./s) 1 / 2 is the r velocity. 
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FIG. 2. Scatterplots of AM reco ii vs AM. The two contours for each mode enclose regions deter- 
mined with signal MC simulation to contain 40% and 85% of signal events, respectively. The signal 
region, defined by (AM, AM recoil ) = (±20 MeV, ±55 MeV), is shown as a box. 



TABLE I. Upper limits on branching fractions of D + and — > /i ± e =F e + at the 90% confidence 
level for a Poisson process [22], where N is the number of D + (or Df) produced in our data, e 
is the signal efficiency, N exp is the number of expected background, iV b s is the number of signal 
candidates, C(A^ bs|-/V exp ) is the 90% confidence coefficient upper limit on the observed events given 
the expected background, and B is the branching fraction or upper limit of the branching fraction 
at 90% confidence level. We increase the upper limits to account for systematic uncertainties by 
decreasing the efficiency, the number of D + (or Df), and the expected number of background 
each by 1 standard deviation. For the D + and — > <fr(e + e~)ir + channels, we have shown both 
branching fractions and upper limits. 



Channel 


N 


e(%) 


^cxp 


iVobs 


C(^obs|iVcxp) 


B 


D+ - 


■> 7r + e + e 


4.76 x 10 6 


33.9 


5.7 


9 


9.3 


< 5.9 x 10~ 6 


D+ - 


■> 7r~e + e + 


4.76 x 10 6 


43.5 


1.3 





2.3 


< 1.1 x 10" 6 


D+ - 


+ K+e + e~ 


4.76 x 10 6 


23.1 


4.9 


2 


3.2 


< 3.0 x 10~ 6 


D+ - 


-> K-e + e + 


4.76 x 10 6 


35.3 


1.2 


3 


5.8 


< 3.5 x 10~ 6 


D+ - 


-> Tr + (ft(e + e~) 


4.76 x 10 6 


46.2 


0.3 


4 


7.9 


(1.7±5;| ±0.1) x 10~ 6 
< 3.7 x 10~ 6 


Df~ 


■> 7r + e + e~ 


1.10 x 10 6 


24.3 


6.7 


6 


5.6 


< 2.2 x 10~ 5 


Df~ 


-> 7r _ e + e + 


1.10 x 10 6 


33.4 


2.2 


4 


6.2 


< 1.8 x 10~ 5 


Df~ 


-> K+e + e~ 


1.10 x 10 6 


17.3 


3.0 


7 


9.3 


< 5.2 x 10~ 5 


Df~ 


-> K-e + e + 


1.10 x 10 6 


27.7 


4.1 


4 


5.0 


< 1.7 x 10~ 5 






1.10 x 10 6 


33.9 


0.7 


3 


6.2 


(0.6±g;| ± 0.1) x 10~ 5 
< 1.8 x 10~ 5 



VI. SYSTEMATIC UNCERTAINTIES 



Possible sources of systematic uncertainty in our measurements are summarized in Ta- 
ble II. Uncertainties associated with upper limits are classified into three categories: uncer- 
tainties due to the normalization (the numbers of D + and Df), the signal efficiency, and 
the number of expected background events. 

Uncertainty in the number of D + (-D+) is estimated by adding contributions from uncer- 
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TABLE II. Summary of systematic uncertainties in D + and Df — > /i ± e T e + decays. Uncertainties 
associated with the branching fraction can be classified as three categories: uncertainties due to 
the normalization (the numbers of D + or Df), the signal efficiency, and the number of background 
events. The columns labeled ir + <fi refer to candidates with <p — > e + e~ decays. 



Source 


7r+e+e 


7T+<}!> 


D+ 
n~e + e + 


K+e+e- 


K~e+e+ 


7r+e+e 


7T + <^ 


D+ 

it e+e+ 


K+e+e- 


K-e+e+ 


Normalization 


2.2% 


2.2% 


2.2% 


2.2% 


2.2% 


5.6% 


5.6% 


5.6% 


5.6% 


5.6% 


Tracking 


0.9% 


0.9% 


0.9% 


1.1% 


1.1% 


0.9% 


0.9% 


0.9% 


1.1% 


1.1% 


PID 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


FSR 


1.0% 


1.0% 


1.0% 


1.0% 


1.0% 


1.0% 


1.0% 


1.0% 


1.0% 


1.0% 


Background suppression 


5.0% 


4.2% 


1.5% 


9.4% 


1.5% 


5.2% 


4.5% 


2.1% 


9.0% 


2.2% 


MC statistics 


0.6% 


0.6% 


0.5% 


0.8% 


0.6% 


0.8% 


0.6% 


0.6% 


1.0% 


0.7% 


Efficiency total 


5.6% 


4.9% 


2.9% 


9.8% 


3.0% 


5.8% 


5.1% 


3.3% 


9.3% 


3.4% 


Number of background 


12% 


68% 


20% 


12% 


25% 


12% 


26% 


16% 


15% 


11% 



tainties in integrated luminosity [12] 1.0% and the production cross section [12] 2.0% (5.5% 
for Df [21]) in quadrature. We assign relative uncertainties of 2.2% to the number of D + 
and of 5.6% to the number of Df. 

There are several sources which can contribute to uncertainty in the signal efficiency esti- 
mation, as listed in Table II. By adding contributions from tracking [12], particle identifica- 
tion (PID) [12, 13], FSR [13, 14], background suppression, and MC statistics in quadrature 
we found total uncertainties in the signal efficiency for each channel range from 3% to 10%. 

We use the number of background events estimated by the MC simulation rather than 
using the sidebands in data. The MC samples, being 5-20 times larger, have higher precision. 
We have evaluated possible systematic bias caused by the use of MC events rather than the 
data sideband by using alternative background shapes, and by comparing the MC predicted 
number to that interpolated from the data sideband. We found no indication of system- 
atic bias; all deviations are adequately explained as statistical fluctuations due to the data 
statistics. We conclude that our MC events reproduce the features of the data backgrounds 
well. We took the statistical uncertainty in the MC simulated number of backgrounds as the 
systematic uncertainty in the expected number of background, as summarized in Table II. 

VII. SUMMARY 

With the complete samples of CLEO-c open-charm data, corresponding to integrated 
luminosities of 818 pb" 1 at E CM = 3.774 GeV containing 2.4 x 10 6 D+D~ pairs and 602 
pb^ 1 at Eqm = 4.170 GeV containing 0.6 x 10 6 D^DJ pairs, we have searched for rare 
(FCNC) and forbidden (LNV) decays of D + and Df mesons of the form h ± e T e + , where 
is either a charged pion or a charged kaon. We found no evidence of signals and set 
upper limits on branching fractions at the 90% confidence level as summarized in Table I. 
Systematic uncertainties in the signal efficiency, the number of D + (or Df) events, and the 
expected number of background events are incorporated by decreasing the numbers used for 
those quantities by 1 standard deviation of the systematic uncertainty on those quantities. 
These results are the most stringent limits on FCNC and LNV for the D + and D+ — > h ± e T e + 
decays to date and the limits in the dielectron channels are comparable to those in the dimuon 
channels [18], but are still a few orders of magnitude larger than the SM expectation [3, 5] 
in FCNC decays. This leaves some room for possible enhancement [2-5] in both FCNC and 
LNV decays induced by non-SM physics. We have separately measured branching fractions 
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of the resonant decays D + — > n + (f) — > 7i + e + e~ and — > n + <f) — > 7i + e + e~ due to their large 
expected contributions to 7r + e + e~ channels. The significance of our measured branching 
fractions is poor at 3.5 standard deviations for D + and 1.8 standard deviations for Df, so 
we have also included upper limits in Table I. Our measured branching fractions of these 
decays are consistent with the products of known world average [18] branching fractions, 
B{D + 07T+ -> e + e-7r+) = B{D + -> 0tt+) x B(<p -> e+e") = [(6.2 ± 0.7) x 10~ 3 ] x [(2.97 ± 
0.04) x 10" 4 ] = (1.8 ±0.2) x 10~ 6 and B(D+ -> <pix + e + e-n + ) = B(D+ 0tt+) x B{(j) -> 
e+e~) = [(4.38 ± 0.35) x 10~ 2 ] x [(2.97 ± 0.04) x 10~ 4 ] = (1.3 ± 0.1) x 10~ 5 . 
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